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Abstract 

A certain vector-tensor (VT) theory is revisited. It was proposed and analyzed as a theory 
of electromagnetism without the standard gauge invariance. Our attention is first focused on a 
detailed variational formulation of the theory, which leads to both a modified Lorentz force and 
the true energy momentum tensor of the vector field. The theory is then applied to cosmology. 
A complete gauge invariant treatment of the scalar perturbations is presented. For appropriate 
gauge invariant variables describing the scalar modes of the vector field (A-modes), it is proved 
that the evolution equations of these modes do not involve the scalar modes appearing in General 
Relativity (GR- modes), which are associated to the metric and the energy momentum tensor of the 
cosmological fluids. However, the A-modes modify the standard gauge invariant equations describ- 
ing the GR-modes. By using the new formalism, the evolution equations of the A-perturbations 
are derived and separately solved and, then, the correction terms -due to the A-perturbations- 
appearing in the evolution equations of the GR-modes are estimated. The evolution of these cor- 
rection terms is studied for an appropriate scale. The relevance of these terms depends on both 
the spectra and the values of the normalization constants involved in extended electromagnetism. 
Further applications of the new formalism will be presented elsewhere. 
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I. INTRODUCTION 



In previous papers, an extended theory of electromagnetism was proposed [1| and devel 



oped 2j-|4|. The evolution of cosmological scalar perturbations was studied in |2|; never- 
theless, the authors assumed that the scalar perturbations of the electromagnetic field do 
not affect metric perturbations evolution. They stated that this assumption holds both in 
the radiation dominated period and in the matter dominated era, that is to say, at any 
moment before dark energy domination. Hence, it affects the choice of the initial conditions 
for numerical integrations, which are always fixed in the radiation dominated era (z ~ 10 8 ), 
when the cosmologically significant scales had superhorizon sizes. Moreover, in paper 2[ , the 
conservation law of standard electromagnetism, V^J M = 0, is assumed, e.g., to get Eq. (2.3) 
from Eq. (2.2). Nevertheless, this law is not an equation of extended electromagnetism, 
whose true conservation law (see below) admits solutions with V^J M ^ 0. This fact is im- 
portant in cosmology, where the law V M J M = implies that only vector modes are involved 
in the expansion of the current J M , whereas the condition V^J^ ^ requires the existence of 
J^ 1 scalar modes. We have introduced one of these modes (see below) in our general calcula- 
tions. We see that, in extended electromagnetism, there are scalar modes in the expansions 
of both the electromagnetic field A 11 (A- modes) and the current J M (J- modes). Since we have 
not convincing arguments proving that, initially, at z ~ 10 8 , all these modes are negligible 
against the small scalar modes of the radiation fluid, we cannot ensure that they do not 
affect metric perturbations. Hence, the A and J scalar modes should not be neglected a 
priori in order to get the metric and fluid initial conditions for numerical integrations. After 
these considerations it seems that the assumption used in j^J is a simplifying condition which 
would require further justification (if it exists). In this situation, it is obvious that a more 
general study of cosmological perturbations is worthwhile. It is performed in this paper, 
where a general complete treatment of the cosmological perturbations is developed in the 
framework of extended electromagnetism. Our approach has various relevant properties: (i) 
it is gauge invariant, (ii) it does not involve approximating conditions, (iii) it involves a 
scalar mode as it is required by the general conservation law of extended electromagnetism, 
and (iv) it uses appropriate scalar modes for the field A 11 which evolve independently of the 
scalar GR- modes (metric and fluid modes). 

This paper is structured as follows: The basic equations of the VT theory are derived 
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-by using variational techniques- in Sec. [TT1 where the energy momentum tensor and the 
Lorentz force are calculated. The theory is applied to cosmology in Sec. IHIt where the linear 
approximation is studied by using standard techniques based on decoupled scalar, vector, 
and tensor modes [§[ {J. An analysis of the equations derived in Sec. II III is performed in 
Sec. IIV[ where the equations satisfied by the scalar perturbations of the vector field are 
numerically solved, and the differences between the equations describing the evolution of 
the scalar modes in General Relativity (GR) and in VT are studied. The Appendix contains 
some estimates in the framework of GR, which are useful in Sec. IIV1 Finally, Sec. IVl contains 
a general discussion and our main conclusions. 

Let us finish this section fixing some notation criteria. Latin (Greek) indexes run from 
1 to 3 (0 to 3). The gravitational constant, the scale factor, the conformal time, and the 
Hubble constant are denoted G, a, r, and Hq, respectively. Whatever function D may be, 
D/ a stands for its partial derivative with respect to the coordinate x a , and Db represents its 
background value. Units are chosen in such a way that the speed of light is c = 1. Spatial 
distances are given in Megaparsecs. 



II. THE THEORY: VARIATIONAL FORMULATION AND BASIC EQUATIONS 

A charged isentropic perfect fluid is considered in the framework of the VT generaliza- 
tion of Einstein-Maxwell theory proposed in 
following action: 

R 1 



1[. The basic equations are derived from the 



16nG 4 



- - A F^F, U + ^V^f + J»A, - p(l 



-g d A x 



where 7 is an arbitrary parameter, R, g^, and g are the scalar curvature, the covariant 
metric components, and the determinant of the g^ v matrix, respectively. The vector field 
of the theory is A^. The symbol V stands for the covariant derivative and we define F^ v = 
V fj,A u — V ' U A^. The electrical current is J M = p q U^, where p q is the density of electrical 
charge and is the four- velocity of the fluid world lines. Finally, for an isentropic perfect 
fluid, one can introduce a conserved energy density p [V At (pt/ ft ) = 0] and an internal energy 
e; so the fluid energy density is p = p(l + e) and the pressure is P = _P 2 (de/dp) (see j^j). 
Some VT theories were proposed in the early seventies (see jf], Q])- All these theories 
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were based on the action: 

/ = (levrGr 1 J (R + uA^R + r J R liV A' i A v - eF^ v + r V^ (1 V# + L m ) 

sFg<£x (2) 

where w, 77, e, and r are arbitrary parameters and L m is the matter Lagrangian, which 
couples matter with the fields of the VT theory. Actions ([I]) and ([2]) are equivalent for 
u) = 0, 2e — 77 = 8irG, t = rj — lQirG'y, and L m = J^A^ — p(l + e). 

According to the variational techniques described in 5] , three fields may be independently 
varied in the action ([!]). These fields are: the vector field of the theory (A^), the flow lines 
of the fluid (U^), and the metric field (57^). 

We first vary the field A^ for fixed flow lines and metric (S A variations). Thus, we easily 
obtain the field equations for the A^ field, whose form is: 

Vi^ = J, + K » ( 3 ) 

where = — 27V M (V • A) with V • A = V M A M . Then, from these field equations one easily 
gets the relation: 

v% = -w; , (4) 

which indicates that the total current J^+J^ is conserved in the theory. This is the conserved 
current associated to the invariance of action ([1]) under the residual gauge transformation 

a; = a^ + with <9^0 = 0. 

In a second step, the flow lines are varied for fixed A^ and g^ u (5 Lr variations) and, 
moreover, the densities p and p q are adjusted to satisfy the equation V At (pf/' t ) = and 
Eq. (j4j), respectively (see j^j). Since the right hand side of Eq. fl4]) does not depends on 
£7 M , the following relation is satisfied 5 u (V l _ l J' J ") = V fJi (5 u J' Jl ) = 0. On account of these 
considerations, the following equations are easily obtained p]: 



{H + P)U^V^U U = -V M P(^ + WIT") + + (V% )A V . (5) 

These equations describe the fluid evolution in the VT theory. The last term is the gener- 
alized Lorentz force, f L , of the theory; hence, we can write: 

fl = V + (V% A )^ (6) 
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Finally, the metric is varied whereas fields and are fixed (5 variations). For this 
kind of variations, Eq (j3J) leads to the relation S g (^—gJ a ) = —5 g (\/—gJ Aa ) and, then, from 
this relation and the identity V^A^V • A)} = (V • A) 2 + A^V^V • A), it follows that the 
Lagrangian densities 7(V M A M ) 2 + J^A^ [involved in action ([T])] and — 7(V A1 A M ) 2 are fully 
equivalent. On account of this fact, 5 g variations lead to: 

GT = 8ttGT^ , (7) 

where is the Einstein tensor and T^ v is the energy momentum tensor of the charged 
fluid plus the electromagnetic field, whose form is: 

T" v = (fM + P)UW + Pg" v + F» a F ua - -g^F^F * 

+2 7 [{A Q V Q (V • A) + ^(V • Af}g^ - A^V"(V ■ A) - A"V^(V ■ A)} . (8) 

This energy momentum tensor is to be compared with that given in taking into account 
differences in the assumed signature. 

Equations ([3D, (EJ), and ((7j) are the field equations of the theory. Any solution of these 
equations satisfies Eq. (j3J), and also the relations V^T^ = 0. These last relations combined 
with all the field equations and Eq. (jSJ) lead to: 

\U p Vp{p) + (/i + P)V P U P ] U a = 2F aP J A p - 2A a V p Jp . (9) 

In the Einstein-Maxwell theory, which is formally obtained in the limit — > 0, Eq ([9]) 
reduces to U P V p(p) + + P)V ' pU^ = (which is identical to Eq. (3.9) in |JJ). 

After this detailed variational study, which has not been previously developed, we are 
concerned with cosmological applications. 

III. BACKGROUND UNIVERSE AND COSMOLOGICAL PERTURBATIONS 

From the field equations of Sec. Dll one easily finds the equations describing a homogeneous 
and isotropic neutral flat universe, in which the line element is 

ds 2 = a 2 [-dr 2 + dr 2 + r 2 d9 2 + r 2 sin 2 9dtf] . (10) 

The vector field, A, and the four-velocity, U, of the cosmic fluid have the covariant com- 
ponents [Aqb{t), 0, 0, 0] and [— a(r), 0, 0, 0], respectively. The energy density is the critical 
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one, whose present value is pso = SHq/8tiG, and the density of charge is Pqsij) = for any 
time. In this background, Eqs. (j3J) reduces to 

S fl = (V ■ A) B = -\[A 0B + 2-Aob] = constant , (11) 
and Eq. (jSJ) leads to the relations 

Pb = ^ = -73|, , (12) 

where quantities p^ and Pg are the background energy density and pressure of the vector 
field A 11 , respectively. The equation of state (fl2j) proves that the energy density of the 
background field A^ play the role of a cosmological constant. In order to have positive 
values of p^ the parameter 7 must be negative. Hereafter units are chosen in such a way 
that c = 8nG = 1; thus, from Eqs. (J7|) and (JS]) one easily get the following basic cosmological 
equation for the background evolution: 

A = a ^PB+P A B) (13) 

a 

~2- + - 2 =a\P B + Pi) (14) 
a cr 

where and Pb are the background energy density and pressure of the cosmological fluid 
(baryons plus dark matter and radiation). Hereafter, w and c 2 s stand for the ratios Pb/ Pb 
and dPs/dps, respectively. 

In next sections, perturbations are described -in the usual way- with the formalism 



summarized in [8j (see also [9|). There are three types of perturbations whose evolution 
is independent during the linear regime. They are the so-called scalar, vector, and tensor 
fluctuations, which may be expanded in terms of the scalar, Q^°\ vector, \ and tensor, 
Q\f , harmonics, respectively. 



A. Tensor perturbations 

There are no tensor modes involved in the expansion of vectors A M and J M ; hence, in 
the VT theory, tensor modes only appear in the expansions of the same quantities as in 
GR (metric and anisotropic part of the stress tensor) and, moreover, they satisfy the same 
equations as in GR. Therefore, we are mainly interested in scalar and vector modes. Metric 
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tensor modes (gravitational waves) and anisotropic stress tensor components evolve as in 
GR, namely, they satisfy the equation: 



(15) 



where li^Qf^ is the tensor part of the anisotropic stress tensor and H@)Qy9 the tensor 
part of the metric (see 8] ) . 

For negligible anisotropic stress, cosmological fluctuations evolving well outside the effec- 
tive horizon (a/a » k, see js, 10]) obey the equation + 2(a/a)H^ ~ 0, whose general 
solution is oc a~ 2 ; hence, is a fast decaying mode. This means that well after 
reheating, e.g. at z = 10 8 , superhorizon scales evolve in such a way that quantity is 
almost independent of time. This fact will be taken into account below. 



B. Vector perturbations 

In the case of a flat background, the vector harmonics can be written as follows jo]: 
Q ± = c 11 exp(ik ■ r), where k is the wavenumber vector. A representation of vectors e + and 



e is 
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ef = (±hk 3 /k - ik 2 )/aV2 , ef = (±k 2 k 3 /k + , ef = ^a/k\/2 , (16) 

where a = (kj + kfj 1 / 2 . 

The first order vector part of the fields and J M may be written as follows: 



(17) 
(18) 



As it is done in standard electromagnetism (Einstein-Maxwell theory), we define the 
covariant components of the electric and magnetic fields (E^ and B^) as follows: 



= F iiU U v , 



(19) 



and 



B, = -S-9Y ll2 ^ vpX F^ 



(20) 
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where quantities e^ vp \ are the Levy-Civita symbols. Then, at first order in perturbation 
theory, one easily finds: 






—aE\ 


-aE 2 


-aE 3 


aEi 





aB 3 


-aB 2 


aE 2 


-aB 3 





aBi 


aE 3 


aB 2 


—aBi 






(21) 



Equations. f[T6"j) - ffT%j) . and ( l2Tj) lead to the following equations in momentum space: 

a a 



(22) 



Hereafter we use the following equivalences: E = (E±, E 2 , E 3 ), B = (Bi, B 2 , B 3 ), and J = 
(Ji, J it J?>)- Moreover, V • X and VAX stand for the ordinary divergence and curl of X, 
respectively. 

By using this notation and Eqs. f[2"2"j) it is easily verified that, up to first order, the 
equation 

a(V A E) + -^-(aB) = . (23) 

is satisfied in position space. Finally, since the ordinary divergence of any vector mode 
vanishes, we can write 

V-B = B ifi = . (24) 
Let us now consider the field equations ([3D, which may be rewritten as follows: 

[In (V=0)] /a F*» + F** /a = JP- 2 7 /°(V ■ A) /a . (25) 

Since V • A and V • J are scalars, they may be expanded in scalar harmonics with no 
contributions from vector modes; hence, the term — 2 , yg l3a ('V ■ A)/ a vanishes in the case of 
vector modes and, consequently, Eqs. ( 1251) reduces to those of Einstein-Maxwell theory. By 
using Eqs. (12 ip and fl25|) . it may be easily verified that -for vector modes and up to first 
order, Einstein-Maxwell field equations and Eqs (|25|) may be written as follows: 

d 



a(VAB) - ^-(aE) = a 2 J 
or 



(26) 

V • E = Et/i = . (27) 

Finally, from Eqs. f TTS]) . (f2"2"|) . and (12"6"|) . the following equation describing the evolution of 
A (1)± is found: 

i(i)± + tf A (i)± = a 4 j(i)± . (28) 



For a given function ./(^(fc, r), the solution of Eq. (128]) gives A^^ik, r) and then, 
quantities E^± and 5 (1)± are fixed by Eqs. (1221) . From these last quantities we may calculate 
E(x,t) and B(x,t) by using the explicit form ( 1T61) of the vector harmonics. The resulting 
E and B quantities satisfy the four Eqs. ( 1231) . ( JiHl) . (1261) . and ( |27j) in position space. 

The equations of this subsection are valid in the VT theory under consideration as well 
as in the standard Einstein Maxwell theory. It is due to the fact that vector modes do not 
contribute to the terms involving 7 either in Eq. (|25|) or in Eq. ([8]), and these terms are 
responsible for all the differences between both theories. The predictions of these theories 
only may be different due to the scalar modes involved in the vector fields and J M , which 
are studied in next subsection. 



C. Scalar perturbations 

For a flat background, the scalar harmonics are plane waves; namely, = exp(ik ■ r). 
The first order scalar contributions to vectors A^ and J M are: 

A, = [a^Q(°\^Qf}, (29) 

J M = [0,a 2 j(°)Qf] , (30) 

where Qf = {-\/k)Qf. 

Since we assume that the universe is neutral up to first order, the component J vanishes 
and, moreover, taking into account that the equation V • J = is not an equation of the 
VT theory, a scalar part a 2 J^Q^ must be included in the expansion of Jj. Equations ( )29|) 
and ( 1301) are absolutely general. 

Equations ( )2T|) and ( )29|) may be combined to get: 

E® = --(ha® + $W) , fl(°) = 0. (31) 

(X 

Similarly, the scalar V • A may be expanded in terms of scalar harmonics; namely, we can 
write: 

V-A=~ B (1 + Z^Q^) . (32) 

In order to calculate = (V • A)^ we must use the relation V ■ A = dA^/dx^ 1 + 
T^ U A V , which involves the Christoffeld symbols. It is then evident that E^ ) depends on the 
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coefficients a^ ^ and (3^ appearing in the expansion of [see Eq. (12 9p ]. and also on the 
coefficients involved in the expansion of the metric components g^, which appear in the 
Christoffeld symbols. If the resulting is used to write Eqs. (1251) . ([7]) and ([8]) up to first 
order in scalar modes, all the equations are coupled among them and their solution only 
have been found under special assumptions (see |2f and Sec. [I]). However, we have developed 
a method, in which we may first solve Eqs. ( |25l) and, then, the solution can be used to solve 
Eqs. ([7]) and ([H]). Let us now describe this method and derive the evolution equations for 
the scalar modes. 

Coefficients and f3^ are not the most suitable ones in order to expand the field 
equations. It is preferable the use of the coefficients and H^, which are gauge invariant 
quantities. In terms of these variables, Eqs. fl25|) reduce to: 

27aS B S<°> = kE^ (33) 

£(o) = _ a 3 j(o) _ 2l ka~ B E^ - ^E® , (34) 

a 

and Eq. (j4]) may be written as follows: 

H(o) + 2 ^(°) + fc 2 »(°) = —a 2 J^ . (35) 

a 2^E B 

Since Eq. (B| is a consequence of Eqs. ( 125]) . which are equivalent to Eqs. (jHJ), Eq. ( 13~5]1 may 

be easily obtained by combining Eqs. ( 133|) and (134|) . Hence, functions H'- -' and .E^ -* may be 

found by solving Eqs. ( 133]) and (I34p for a given J 1 - ) plus initial values of 5^ and E^ and, 

then, the resulting and E^ - 1 functions and the chosen will satisfy Eq. ( 135]) . 

Equation ( 134"]) may be easily derived from Eqs. ( 133]) and (135]) and, consequently, we may 

also proceed as follows: Eq. (135 p is solved for a given and initial values of S(°) and 

and, then, the resulting solution S^ ^ is used to get function E^ by using Eq. ( 133]) . 

Obviously, Eq. ( 134"P is satisfied by the H^ ) and E^ functions we have found by solving 

Eqs. 035]) and 033]). 

It is worthwhile to point out that Eqs. ffT5]) and ( 135]) have the same form. In fact, if we 
replace H^p by and IiP by —kJ^/2Pg^f in Eq. ( 1151) . the resulting equation is identical 
to Eq. ( 133]) . Condition J(°) = in Eq. ( 133]) is equivalent to condition Wp = in (ITS]) . 
Hence, some previous conclusions about the evolution of gravitational wave modes would 
be also valid for the E^ evolution; in particular, for j(°> = and superhorizon scales, the 
relation E^ ~ holds. 
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In extended electromagnetism, function does not vanish a priori. Condition J(°) = 
implies the relation V^J^ = in position space, but this relation is not a basic equation of 
the theory. By this reason, function J^> is included in the equations derived in this section; 
nevertheless, there is no -by the moment- any physically motivated rule to build up this 
function. Anyway, for a given (including the cosmological possibility = 0), Eqs. (|33|) 
and ( l34l) involve: the constant 7 < 0, the wavenumber k, the scale factor, functions S^ * 1 and 
E(°\ and their first order derivatives. These equations do not involve scalar perturbations 
associated to the metric and the energy momentum tensor. They may be easily solved for 
given values of 7 and k, and initial values of and (alternatively we may solve Eq. f l35|) 
for initial values of H^ ) and H^ )). The resulting function H^ ) appears in the expansion of 
Eqs. ([7]) and (jHJ) in scalar harmonics (see below). 

The gauge invariant formalism described in [8[ is used in this paper; namely, the metric, 
the four- velocity, and the part of the energy momentum tensor (jSJ) being independent of the 
parameter 7 are all expanded as follows (in the flat case): 

g 00 = -a 2 (l + 2AQ ( -% 9oi = -a 2 Qf\ 
9ij = a 2 [(l + 2H L Q<V)6 ij + 2H^Q^] 
U t = av^Q?\ p = p B (l + 5Q^) 

Tij = P B {1 + KLQ^Sij + Ps^Qlf ■ (36) 

Any other quantity as, e.g., Uq, Toi, and so on, may be easily written in terms of the 
coefficients involved in these equations (see 8j), which may be combined to build up the 



following gauge invariant variables: 



rC 



77 = (WTT L — C s 5)/w , V 

k a 

k k a k z V a 



k k a k 2 \ a 

1 „m 1 a 1 a 

3 1 k a k 2 a 



1 ~ cm . 1 a ~ f ri\ T v(0) ' 

r 

^ = ^ + ^4 0) + 7-^ (0) -A-4 0) (37) 



The complementary part of the energy momentum tensor (jSJ); namely, the part depending 
on 7 may be easily expanded in terms of scalar harmonics. The resulting expansion involves 
the variable S^ ^ and its first order time derivative (or equivalently and E^). 
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In order to expand Eqs. (|7j) and (JSJ) and the relation V y T^ LU = 0, we use the same gauge 
invariant potentials and variables as in jg]) (see above), plus the gauge invariant variables 
S(°), E(°\ and The resulting equations reads as follows: 

^OB^(o)" 



2/c 2 
a 2 



®H = pB^r 



2 7 : 



■B 



3— Ann + 



a- 1 



+ 



k 2 
a 2 

v(o) + d L v (o) = k $ A+ k ,J 



PbTiP 



c s e m + wrj) 



2wk 



-n 



(0) 
T 5 



(38) 

(39) 
(40) 



a 1 + w 3(1 + it)) 

(p B a 3 e m y= -ka\p B +P B )v^-2a 2 aP B U^-ka 3 A 0B J^-3 1 a 3 E B A QB (p B +P B )E^ . (41) 

Eqs. (1551) to (T4T1) plus (JT5J and (fT4"j) may be combined to get the following equation for the 
evolution of e m : 



(pBa 3 e m )"+ (1 + 3c 2 )-(p B a 3 e m ) 



k2c l - ^(Pb + ^Pb)^ 2 (Pi?a 3 e m ) 



-k 2 (P B a 3 r)) - 2d(P B a 2 4 0) )-+ ^ 2 (P B a 3 4 0) ) 



+2p B a 2 



w 



l + < - 



:p B a 3 4 0) ; 



-2 7 a 3 ~ B (p B + P B )(2A 0B S(°) - S B a 2 S(°)) 



ka 3 A 0B j^ + 



(2 + 3c 2 s )-A 0B - a 2 E B 



J®) . 



(42) 



For EP' = 0, J^ ^ = 0, and p B = 0, Eqs. (138|) to ( 142|) reduce to the equations derived by 
Bardeen in the fiat case |s|. For E^ = 0, J( 0) = 0, and p B = p A ^ 0, Eqs. (1381) to (1421) 
describe fluctuation evolution in a standard flat universe with a cosmological constant whose 
energy density is p\. Finally, if p B = p\, and the two functions and E^ do not vanish 
at the same time, Eqs. ( 138|) . ( 14T|) . and (142]) contain new terms, which modify the equations 
describing perturbation evolution in flat universes with cosmological constant. 



IV. ANALYZING THE BASIC DIFFERENTIAL EQUATIONS OF EXTENDED 
ELECTROMAGNETISM 

In this section, it is assumed that the condition = holds in cosmology, which is 
equivalent to assume the well known conservation law of Einstein-Maxwell theory (V M J M = 
0). Under this arbitrary assumption, we study the background equations (irTT) -(lT4l. and 
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the equations ( 133|) - (!35|) . and (138|) - (|42|) describing the evolution of the first order scalar 
perturbations in momentum space. 

We begin with the background equations. Function p B (r) and P B {j) are given by the 
formulas 

pB = PBvo{1 + Zf + PBm0 (l + zf , P B = PBr0 (l + zY/3 , (43) 

where p Br0 = 8 x 10~ 34 gr/cm 3 and p Bm o = 0.2726p c are the present energy density of 
radiation and matter 

Moreover, the baryon density is assumed to be p B bo = 0.0461p c , and the value of the 
Hubble constant is Hq = lOOh Kms~ 1 Mpc~ 1 with h = 0.704. All these values are compatible 
with a certain version of the concordance model (see [3]). In this model, the dark energy 
density, p B , is easily obtained from the relation p Bm o + PBro + P B = which is valid in 
flat backgrounds. 

Equation ffl3|) governing the evolution of the scale factor may be numerically solved for the 
above parameters; thus, the evolution of the scale factor is obtained. The resulting function 
cj(t) is necessary to study the remaining equations of the theory; namely, Eqs. fl33|) - fl35|) 
and also Eqs. (|38|)-(|42|). 

From Eq. ffT2l one easily gets the relation E B = ±|7|~ 1 / 2 (p^) 1 / 2 (negative 7). Once the 
value of p B = pa is fixed (see above), this last relation leads to E B oc S^^l^ 1 / 2 , where S gn 
only may take on the values +1 and —1. For a given value of I7I, only the absolute value of 
E B may be obtained (its sign is arbitrary). On account of Eq. fTlTT) . we may also write the 
relation A 0B oc S^^l -1 / 2 . In the background, quantities I7I and S gn remain arbitrary. 

Hereafter, D m stands for the initial value of quantity D whatever it may be. Let us now 
consider Eqs. (l33l -(|35l ) . We first solve Eq. ( 1351) by using initial values S(°) m and at 
redshift z = 10 8 . At this high redshift, the cosmo logical scales of interest (see below) are 
superhorizon ones; hence, taking into account the similarity between Eqs. (1151) and fl35l) and 
the comments in the last paragraph of Sec. IHI A\ the condition = assumed in this 
section allows us to take = 0. Only the initial value of may be appropriately chosen 
to integrate Eq. ( 1351) . The values of I7I and S gn are fully irrelevant to perform this integration. 
Hereafter, numerical calculations are performed for the spatial scale L = 3 x lO 3 ^" 1 Mpc, 
which reenters the effective horizon at present time (tq). Its wavenumber is k ~ 1.47 x 10~ 3 . 
This scale is useful for normalization in GR (see the Appendix). Equation (1351) has been 
solved for the wavenumber k with the initial condition H^ m (L) = 10~ 4 in position space. 
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In order to write the corresponding initial condition in momentum space, we use the well 
known relation jisj 

(\X(x)\ 2 ) L — k 3 (|X(A;)| 2 ) /2tt 2 , (44) 

where X is an arbitrary quantity. Evidently, this relation must be particularized for X = 
«;(°) m ; L — L, and k = k to calculate S(°) m (/c). Either this last initial quantity (in momentum 
space) or E^ m (L) (in position space) may be seen as a normalization constant. The S^ ) 
spectrum would be necessary to derive the initial value of S^°^(fc) for k ^ k. 

The solution of Eq. fl35|) plus Eq. f )33|) allow us to calculate function which is pro- 
portional to Sgnl'y] 1 / 2 . 

Numerical integrations have given the functions A 0B (z), E^(k,z), and E^(k, z) rep- 
resented in Fig. [TJ These functions correspond to |7| = 1, S gn = +1, k = k, and 
«;(°) m (Z) = 10~ 4 . Since the dependence of these functions in terms of the parameters 
I7I, S gn , and 5(°) m (fc) is known (see above), Fig. [1] contains complete information about 
the scalar modes associated to the field A 11 for the scale k. The same may be done for any 
linear spatial scale with the help of an appropriate spectrum for S(°) m . Let us now study the 
Eqs. (I38l) - (14"2"|) describing the evolution -in the framework of extended electromagnetism- 
of the scalar modes appearing in standard GR cosmology. Equations (I38I) . (14T1) . and (142]) 
contain the terms: 

6(r, k) = -2724(3-%, + 2 B )2(°) + SW] , (45) 

. \ CI, / (X 

6(r, k) = -3 ia 3 E B A 0B (p B + P fl )H(°) , (46) 

6(r, k) = -2 ia s E B (p B + P B ){2A QB m - ~ B a 2 ~W) , (47) 

respectively. These terms -which appear in extended electromagnetism but not in Einstein 
theory with cosmological constant- may be calculated, for the scale k, by using the integra- 
tion data used to build up Fig. CD and, then, these terms may be compared with appropriate 
terms involved in GR equations (for the same wavenumber). 

Taking into account that S B and A 0B are proportional to S gn \^\~ l l 2 and also that E^ ) 
does not depend on S gn and I7I, it is trivially proved that quantities £1, £ 2 , and £3 are also 
independent of S gn and |7|. For the wavenumber k, these quantities are proportional to the 
number E^ m (L). 
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As it follows from Eqs. (158]) . ( )4T|) (142]) . quantities £1, £2, and £3 are to be compared with 
the GR values of the terms 

T 1 {T,k) = p B e m , (48) 

T 2 (r,k) = -ka 3 (p B + P B )v^ , (49) 

T 3 (r, fc) = [^(Pb + ^)« 2 - k 2 c^ (p B a 3 e m ) , (50) 

respectively. After the estimation of T 1; T 2 , and T 3 in standard cosmology (based on GR), 
the three functions r^r, k) = |£i(r, fc)/Tj(r, may be calculated. Evidently, for very small 
r, values, GR and extended electromagnetism would lead to the same differential equations 
for the evolution of the GR scalar perturbations, whereas values of the order of 10~ 3 or 
greater would suggest relevant differences with respect to GR. If differences are expected, 
the matter power spectrum P(k) and the angular power spectra of the CMB should be 
accurately estimated by using numerical codes as CMBFAST 17] and CAMB [14]. These 
accurate calculations -under general enough initial conditions- are beyond the scope of this 
paper. 

The Ti ratios will be estimated for the wavenumber k. The corresponding spatial scale is 
useful for normalization in GR, which is necessary to estimate the functions Ti(r,k). The 
method used for the estimation of these functions and for normalization in standard GR 
cosmology are described in the Appendix. 

The three functions J"j(z, k) are represented in Fig. [2] for ES°' 1,n (L) = 10 -4 . Two panels 
(left and right) show the evolution of each ratio ?y The evolutions of the three ratios are 
similar. From z = 10 8 to z ~ 10 2 , the chosen spatial scale is well outside the effective horizon 
and all the ratios increase without oscillations (see left panels); however, for z < 10 2 , there 
are oscillations whose amplitudes grow as z decreases (see right panels). It is due to the fact 
that our spatial scale and that of the effective horizon come near as the redshift decreases 
(they have been chosen to be identical at z — 0). The maximum values of the ratios are 
n ~ 6.3 x 10~ 4 , T2 — 2.8 x 10 -4 , and r 3 ~ 2.6 x 10~ 3 . These maximum values are all reached 
close to z — 0. They are small enough to ensure that, for S^°' )m (L) = 10~ 4 , the chosen scale 
evolves as in GR. 
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V. DISCUSSION AND CONCLUSIONS 



An exhaustive variational formulation of extended electromagnetism has been presented 
in Sec. [Ill In particular, the energy momentum tensor of the vector field A^ 1 has been 
found. This tensor has been analyzed to conclude that the last term of Eq. flS} has not 
the same sign as in previous calculations Our sign appears as a result of the relation 
8 g {y/—gJ a ) = — 5 g (\/^gJ a ), which is consistent with the conservation law of J^ + Jf}. This 
sign is only compatible with a negative constant 7 (see Sec JIIIj) . Our formulation also leads 
to Eqs. ([6]) giving the components of the Lorentz force and to Eq. (Q. 

We have emphasized that, for vector and tensor linear perturbation of the Friedman- 



ly equivalent to Einstein- 
, our detailed but brief 



Robertson- Walker universe, extended electromagnetism is fu 
Maxwell theory. Although this fact was already known (see 
study of vector and tensor perturbations is necessary to describe our approach to the evolu- 
tion of scalar modes and, moreover, this study exhibits novel aspects. Note that: (i) our brief 
considerations about tensor perturbations have been used to study the evolution of some 
scalar modes [analogy between Eqs. (|15p and (j3"5|) ]. (ii) our analysis of vector perturbations 
leads to new relations in momentum space derived from Eqs. f )16p and, (iii) this analysis 
suggests the definition of the gauge invariant quantity E^°\ which plays a crucial role in our 
description of the scalar modes. 

A new approach to deal with the evolution of the scalar modes in extended electromag- 
netism has been described. It is gauge invariant and fully general. This formalism leads to 
Eqs (I33|) - (I35|) and Eqs. (I3"51) - (H2I) . which are the most important findings of this paper. 

The scalar modes associated to the vector field are assumed to be 5^ and E^°\ 
Then, in Sec. IIII C\ it is proved that the evolution equations of these gauge invariant modes 
[Eqs (!33|) -( l35|) ] do not involve any other scalar mode (excepting J^). These simple equations 
may be easily solved by using standard numerical methods for any given function 

In Sec. IIVI where the simplifying condition = is assumed, the solution of Eqs (133jl - 
(I35p may be found by using the initial value of H*- -* at z ~ 10 8 , whereas the initial values 
of both S(°) and E° [related by Eq. 032])] may be neglected. The solution corresponding to 
a particular wavenumber k (reentering the effective horizon at present time) is presented in 
Fig. [TJ There is no problem to integrate Eqs ( |33]) -( |35]) for any other scale. Since we have the 
numerical solutions of these equations, the correction terms & appearing in Eqs. ( )38l) -( l42l) 
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may be numerically treated as known functions of k and r and, consequently, these last 
equations only involve -as unknown functions- the scalar modes appearing in GR. Their 
numerical solution has not been found in this paper; where the correction terms £j have 
been compared with appropriate terms of the GR equations (Tj) for suitable values of k 
and S(°) m . For the chosen scale k and S(°) m (L) < 10~ 4 , the scalar modes would evolve as in 
GR. The spectrum of and the standard power spectrum P(k), at initial time, would be 
necessary to perform similar comparisons for all the cosmological scales. If all the scales are 
found to evolve as in GR, both theories are equivalent from the cosmological point of view; 
on the contrary, some appropriate code, as e.g., CMBFAST or CAMB, may be modified to 
estimate -in the VT- the angular power spectrum of the CMB, the matter power spectrum, 
and so on. 

In the background, functions Aqb{t) and H#(t) cannot be fully fixed. Both functions 
are proportional to S^^l -1 ' 2 , but these parameters are arbitrary. Moreover, function 5^ 
and the correction terms £j defined in Eqs. (j45p -(]4"T ]) are independent of parameters S gn and 
|7|. Hence, cosmological considerations cannot fix the values of these parameters. It is not 
surprising, since it is well known (see Q) that the theories based on the Lagrangian (TSJ) with 
J M = (no currents) cannot completely fix the vector field. 

Condition = has been assumed to be valid in cosmology; nevertheless, this condition 
is not strictly required by extended electromagnetism. The question is: What would be a 
scalar j( )-current in cosmology? More research about this scalar mode and its meaning is 
being carried out. 

Finally, let us discuss in detail the fact that our energy momentum tensor and that 
found in 2( have opposite signs. The possible consequences of this difference deserve special 
attention. 

If the relation 7 = 2£ is satisfied, our Lagrangian (with 7) is identical to that used in Q| 
(including £). In spite of this fact, opposite signs appear in the energy momentum tensors. 
As it is explained in Sec. [TT] and summarized in the first paragraph of this section, our sign 
is obtained -from the common Lagrangian- with right variational calculations based on the 
true conservation law of the theory. Equation (4) is actually the conservation law satisfied in 
extended electromagnetism; however, if we take W M = 0, which is not an equation of the 
theory, but the conservation law of standard electromagnetism, the opposite sign is easily 
found in the resulting energy momentum tensor. This sign is not right. 
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For an arbitrary positive £ value and for the corresponding negative value 7 = — £/2, 
the energy-momentum tensor in is identical to our energy-momentum tensor. Hence, the 
Einstein equations are also indistinguishable for these values of £ (positive) and 7 (negative). 
However, the field equations of the vector field are different for the same values; namely, 
for 7 = -£/2. 

Equations may be written in the form VF^ = J M — 27V /J (V-A), where V-A = V^A^, 
and these field equations are to be compared with the equations V v F tlv = J M — £V M (V ■ A) 
appearing in Qj. This comparison shows that the terms — 27V M (V • A) and — £V M (V • A) 
-which modify the field equations of standard electromagnetis- have opposite signs for 7 = 
— £/2. Since the resulting A M field equations are different, distinct predictions seem to be 
unavoidable in general nonlinear applications of extended electromagnetism. The discussion 
of these nonlinear cases is beyong the scope of this paper, where we are concerned with 
cosmological linear applications of the theory. 

Actually, both signs lead to the same conclusions for first order perturbations of 
Minkowski and Robertson- Walker space-times, namely, for the cases considered, e.g., in 
1^ | and 20]. It is due to the fact that, according to Eqs. (TTT1) . ( Tl2|) . and ( 1351) . the signs of 
"Eb and = (V ■ A)^ /E# are arbitrary and, consequently, for 7 = — £/2, these signs may 
be chosen to make identical the scalar parts of the terms — 27V M (V • A) and — £V^(V • A). 
Thus, the linearized A^ field equations derived in our paper become equivalent to those 

n n n 

of previous papers [2J, H9|, 12Q[. Since the energy momentum tensors are also identical for 
7 = — £/2, the modes of positive energy coincide, and the conclusions of papers [3] and [3] 
(with £ > 0) may be also obtained here for 7 < 0. However, only our signs are right and, in 
general, only our A^ field equations should be applied in nonlinear cases. 
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Appendix: Estimating the T, functions in GR 



Since WMAP observations strongly suggest that cosmological perturbations are adiabatic, 
only the case rj = is considered in this section. 

For cosmological perturbations evolving outside the effective horizon (k < aH/2-ir), the 
evolution is essentially independent of the microphysics. This means that the anisotropic 
stress due to neutrinos may be neglected (i ffi = 0); an d also that, in spite of the tight 
coupling between photons and baryons (see 10j) at z > 1100, the transfer of energy and 
momentum between these two species may be forgotten and, consequently, the corresponding 
fluids may be treated as independent. This means that, for super horizon scales, a good 
enough estimate of functions Tj may be done by solving Eqs. (|3"8]) - P3]) for IT^/ = 77 = 0. In 
this case, Eqs. (j3g)) -(P ]) lead to: 



^ m + (l + 3c 2 s )-4/ m + 

CL 



k 2 c 2 s -Up B + P B )a 2 ^ m = , (A.l) 



2 

y m = -ka 3 (p B + P B )v^ , (A.2 
kc 2 1 



y(0) + ^(0) = I 

a a 



(p B + P B )a 2 2k 



*m , (A.3) 



= 2^ • (A ' 4) 

where \l/ m = p B a 3 e m . This system of equations must be solved together with the background 
equations for appropriate initial conditions at z = 10 8 . 

The background is a flat universe with cosmological constant. The energy densities 
of matter, radiation and vacuum correspond to the concordance model (see above). The 
background differential equations may be easily integrated to get a(r), p B (r), and P B {r). 

The integration of Eqs. fjA. lj) — (lA.3j) only requires e^J and t>i°' >m at z = 10 8 . In fact, from 
the first of these values one easily obtains and the initial value of \l/ m may be then 
obtained by substituting vi°^ m into Eq. (IA.3I) . The second order differential equation flA.lj) 
may be integrated by using and 4f™. Function ^ m (r) is then known and vi°^ m may be 



T 



used to solve Eq. flA.3j) and get function ufi 

Since e m and vi are gauge invariant quantities, their initial values may be calculated in 
any gauge. We have used the synchronous gauge to perform this calculation. For super- 
horizon scales, equation (96) of reference [l(| may be used to easily get the following initial 
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conditions: 



S in 



,(0)m 



'-C(kr 



in\2 



(0)m 



0, v. 



V 



(0)in 



H 



3 

~Ck 3 (r in ) 3 , 
-6C(1 + 1)(^) 2 , 



(A.5) 



where the conformal time r m is that corresponding to the chosen initial redshift z = 10 8 , 
C is a normalization constant, and the subscripts 7, b, and c stand for photons, baryons, 
and cold dark matter, respectively. The fluid formed by these three components has the 
following density contrast and peculiar velocity [lo| : 



5 = (pBbSb + PbcSc + Pb^Mp 



B 



,(0) 



[{P3b + P Bb )vl 0) + (pbc + PbcH 0) + (p Bl + P Bl )v^)}/{p B + P B ) ■ (A.6) 



,(0) 



(0)^ 

7 



The initial values of e m and may be easily calculated taking into account Eqs. (131 
(]A.5[) . and (IA.6[) . Equations. (lA.lj) - flA.3j) may be then solved. 

The estimation of quantities Tj requires normalization. The question is: How can we find 
a good enough value of the normalization constant C? It is well known that, for superhorizon 
scales, the gauge invariant quantity 

1 



2 $ g + (aHy^H 

c = 3 TT^ + ® H 



9\aHJ 1 



w 



is conserved 



15] and, moreover, at horizon crossing, the relation 

c 



S(k,r) = 0(l): 



w 



(A.7) 



(A.8) 



is satisfied, where 0(1) is a number of order unity (see 15j and references cited therein). 
Then, normalization may be achieved as follows: in a first step, the matter power spectrum 
at present time P(/c,tq) is obtained, by using CMBFAST, for a certain version of the con- 
cordance model (see [16| for details). The resulting spectrum is represented in Fig. [3j From 
it, we can estimate P(k,To). In a second step, Eqs. ()A.lj) - flA.3j) are numerically solved for 
the scale k and for an arbitrary C value and, then, by combining Eqs. (IA.2j) . (IA.4[) . and 
( 1A.7I) . the function ((k, r) may be easily calculated. Finally, in a last step, the value of the 
normalization constant C is fixed. It is done by using Eq. (1A.8jl to calculate 8(k,To), and 
taken into account that the resulting 5(k, t ) quantity must be identical to P lj/2 (k,T ) for 
the right C value. 
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The 5(k, To) value obtained from the spectrum of Fig. [3] and Eq. (144jl may be easily used 
to estimate the contrast S(L, r ) in position space, the resulting value is close to 10~ 3 as it 
is expected for this scale. 
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FIG. 1: Functions A 0B (z) (top), Z( \k,z) (middle) and E^(k,z) x 10 8 (bottom) in terms of 
log(l + z), for L = 3000/1" 1 Mpc, | 7 | = 1, S gn = +1, and E^ in (L) = 10" 4 . 
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FIG. 2: Left: functions log[ri (k, z)\ (top), log[r2(fe, z)] (middle) and log[r3(fc, z)] (bottom) in terms 
of log(l + z), from z = 10 8 to z = 0. The wavenumber k is the same as in Fig. [TJ Right: functions 
r±(k,z) x 10 4 (top), r2(k,z) x 10 4 (middle) and r^(k,z) x 10 (bottom) in terms of log(l + z), from 
z = 10 2 to z = 0, for the same wavenumber as in the left panels, 
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FIG. 3: matter power spectrum -estimated with CMBFAST- for the chosen version of the concor- 
dance model. 
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